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Abstract 

We study the holographic flow of anomalous conductivities induced by gauge and gravita- 
tional Chern-Simons terms. We find that the contribution from the gauge Chern-Simons 
term gives rise to a flow that can be interpreted in terms of an effective, cutoff dependent 
chemical potential. In contrast the contribution of the gauge-gravitational Chern-Simons 
term is just the temperature squared and does not flow. 
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1 Introduction 



Anomalies appear in the context of relativistic quantum field theories. In four dimensions chi- 
ral anomalies [1] involve triangle diagrams with either only vector currents or vector currents 
and the energy momentum tensor, in which case one speaks of a (mixed gauge-) gravitational 
anomaly [2]. They are responsible for the breakdown of a classical symmetry due to quantum 
effects. If the symmetry is local anomalies impose severe restrictions on the structure and 
definition of gauge theories (for comprehensive reviews on anomalies see [3j). In the case of 
a symmetry generated by Ta, and considering only right-handed fermions, the presence of 
a chiral anomaly in vacuum is encoded in a non-vanishing cIabc — |tr (Ta{Tb, 7b}). The 
corresponding parameter in the case of the gravitational anomaly is bA = tr (Ta)- 

Recently, it has been pointed out that at finite temperature and density, anomalies are 
responsible for the appearance of new non-dissipative transport phenomena [U EJ E] • In the 
chiral magnetic effect an external magnetic field induces a current parallel to it 



J M = o B B» (1) 

where <jb is the chiral magnetic conductivity and B p = ^e tJ-upx u u F p x- A second effect is 
the chiral vortical effect. It refers to the creation of a current parallel to vortices in the fluid 



^ = avU » (2) 

with u p = e^ vpX u v d p u\ being the vorticity vector and the fluid four-velocity. 

The contribution of the gravitational anomaly to these transport coefficients was first 
obtained in a weakly coupled gas of chiral fermions in [TJ. A holographic model that confirmed 
these findings at strongly coupling was developed and studied in [8]. 

The chiral magnetic and the chiral vortical conductivities can be calculated from first 
principles via the Kubo formulae [9] (Latin letters denote purely spatial indices) 

a B = lim -L Y^e abc {j a J b ) (u = 0,p) (3) 

a,b 

av = i% i £ €abc < rT ^ {u = °' & (4) 

a,b 

There are related transport coefficients for the energy current T 0fl 

T° M = a e B B p , (5) 
T° M = a e v u p . (6) 
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They are calculated via the Kubo formulae 



(1.6 



(7) 




a.b 
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The role played by the gravitational anomaly was further studied in [10] , where an ideal 
Weyl gas in arbitrary even dimensions was considered. This leads to a generalization of the 
anomalous conductivities valid for any (even) dimension and an expression that relates the 
anomaly induced transport coefficients to the anomaly polynomial of the Ideal Weyl gas. 
Furthermore, in [TTJ, a definition for the local entropy current for higher-curvature gravita- 
tional theories was proposed and the Fluid/Gravity correspondence was applied to compute 
the first order conductivities in the presence of the gravitational anomaly. 

Within the gauge-gravity duality the running with the holographic coordinate can be 
interpreted as a type of renormalization group (RG) flow in the dual field theory [12]. The 
first application of this holographic flow to transport coefficients is [13] where it was shown 
that the electric conductivity and the shear viscosity have a trivial flow. It is known now 
that some of the transport coefficients indeed present a non-trivial flow (see, for instance 
|14j). The extension to finite chemical potential has been studied in j!5j . 
Recently, there is a renewed interest in this subject due to the explicit holographic construc- 
tion of the Wilsonian Renormalization Group [16], which has made possible to show that 
multi-trace deformations in the effective action are induced after integrating out high energy 
modes. Finally, in [T7] it was pointed out that all the apparently different frameworks used 
over time to study the holographic flow are actually equivalent. 

It is natural to analyze the holographic flow of the anomalous conductivities as well. In 
this paper we present several approaches to compute the different flows and show that all 
the methods lead to the same results as expected. In Section 2 the setup is presented. We 
interpret the holographic flow, defined as in [IB], as a cutoff flow that arises by varying the 
holographic cutoff at finite holographic coordinate value r = A. Section 3 studies the case of 
gauge fields without taking into account the backreaction onto metric perturbations. Then 
we generalize the previous approach to include also the metric perturbations and present 
a non-covariant method to calculate the flow equations for the retarded Green's functions. 
Somewhat surprisingly we do find a non-trivial flow but give it a natural interpretation as 
a cutoff flow of an effective chemical potential. The flow of the correlators is computed by 
explicitly solving the equations of motion for a system restricted to live between the black 
brane horizon and a hyper surface placed at finite r = A, which acts as a boundary. The 
section ends with a discussion regarding the compatibility of the results so obtained with the 
flow equations. It is shown that both approaches are in agreement. In Section 4 the attention 
is focused on the gravitational anomaly. We discuss subtleties concerning the definition of 
a Dirichlet problem and the necessity for the inclusion of a boundary term to ensure that 
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the correct form for the divergence of the current is found. Contrary to what happens for 
the chemical potential, we find that the temperature term stemming from the gravitational 
anomaly does not flow. 

We conclude in Section 5 with a discussion. 

2 Setup 

Lets show how transport coefficients flow with a variation of the holographic cutoff scale. 
We define the theory with a cutoff as: 




We consider this theory in a general black brane background of the form 



ds 2 = -g tt dt 2 + g rr dr 2 + gad{x % ) 2 . (10) 

We assume that the above metric has an event horizon at r = r# and that every com- 
ponent depends only on r. The boundary is placed at r = A. The metric is also assumed 
to be regular except at the horizon and possibly in the limit A — > oo. The current of the 
holographic dual field theory is 



r=A 



In the gauge A T 
gauge field reads 



the x-component of its variation due to a small perturbation of the 



(12) 

r=A 

where a = da/dr is the r-derivative of the aforementioned perturbation. 
We define a(x,r) as a(x,r) = ^Xa^(x), so that it is normalized at the boundary to 
a(x, A) = a^(x) and a(r) solves the radial wave equation 



J x = -l-V^gg xx g rr a{x,r) 



o(r) + \a{r) (g tt g tt + f*g„ - g rr g rr ) + g rr {u 2 g tt - k 2 g xx ) a(r) = (13) 

On the other hand, we define the electric conductivity at the boundary as J x = <je(A.)E(A), 
where E(A) = —icoa^ is the external applied electric field. Comparing this to equation (|T2|) 
we conclude 
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<te(A) 



e 2 u 



-99 xx 9 rr 



air] 



a(A) 



(14) 



Varying the cutoff A — > A + dA we find the for the differential of the electric conductivity 



(2cte(A) — z 



e 2 o; 



d 
dr 



-gg xx g rr 



da(r)/dr 
a(r) 



(15) 



This equation shows that we can study the flow of the transport coefficients with the cutoff 
reformulating it as the evolution with respect to the coordinate r, by formally identifying r 
with A. 

We can use now the equation of motion for the perturbation a(r) and the definition of 
the conductivity (TBI) to derive the flow equation 



daE^uj, k) 
dA 



-ioj 



-g 



~grrgxx&E 



~Ag xx 



k' 2 

g"" + —g" 

UJ 2 



(16) 



This the flow equation first derived in [13]. It can be solved by demanding infalling 
boundary conditions on the horizon. In particular the flow for the DC conductivity turns 
out to be trivial &e = 0. In this case the electric conductivity is completely determined by 
its value on the horizon via the membrane paradigm 



ere (A) = a E (r H ) 



(17) 



3 Flow of anomalous conductivities 

We will apply now the strategy outlined before to the anomalous transport coefficients. Two 
models will be considered. First we discuss a model in which we neglect the backreaction 
of the gauge field fluctuations on the metric. We will study the interplay between two U(l) 
symmetries which we call vector and axial ones. This allows to model the chiral magnetic 
and the chiral separation effect. A second model will use only one anomalous U(l) symmetry 
but we will also include the backreaction onto the metric. This allows to study also the flow 
of the chiral vortical conductivity and the flow of the anomalous transport coefficients related 
to the energy current. 



3.1 Vector and Axial symmetries 

We will apply the aforementioned strategy to the chiral magnetic conductivity [18] . Its 
proper definition requires the interplay between a vector like U(l) symmetry and an axial 
U(l) symmetry. Holographic models have been investigated in [TjJl 1201 EI]- The model 
allows for the definition of the chiral magnetic conductivity and axial conductivities involving 
external axial magnetic fields. Its action is given by [21] 
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rpV ttiMN _ r?A rpMN _, 

<-> — I V y\ /\ 9 MN V a „2 MN A "r 



4# v 4<£ 

+ - e ^PQ^ M + 3F V pF&) } (lg) 

where V stands for 'vector' and A for 'axial'. The Lagrangian contains two Maxwell 
actions for vector and axial gauge fields and a particular choice of Chern-Simons term. In 
what follows, we will stick to the notation of [2T]; concretely, we define the epsilon symbol 
as e(ABCDE) = —y/—ge ABCDE , with e(rtxyz) = 1 (r corresponds to the fifth coordinate). 

From the boundary term of this action, after perturbing both the axial and the vector 
gauge fields, we obtain an expression for the boundary theory currents 



(19) 

r=A 

(20) 

r=A 



where e pupX = e rfJ-upX . The coefficients in front of the Chern-Simons terms are crucial to 
ensure that the vector current is non- anomalous D^Jy = 0. The axial current is anomalous 
D^Ji? = -| e ^ A (SF^F^ + F£,Ffc) [21]. Comparing with the standard result from the 
one loop triangle calculation, we find k = for a dual strongly coupled SU(N C ) gauge 

theory for a mass less Dirac fermion in the fundamental representation. Note also that both 
currents are invariant under vector gauge transformations but not under axial gauge trans- 
formations. 

The equations of motion for the gauge fields are 



jV N F a +- 



+ ^e** Np Q R (F£ P F£ R + F^ p Fq R ) = (21) 



\V„F$ M + 3ne MNP( l R (F£ P FV R ) = (22) 
9v 

In order to study the flow of the conductivities with the fifth coordinate, we will proceed 
as follows: 

• We introduce an axial and vector perturbation of the gauge fields 

A M = A < $ + a M (y,t,r) (23) 
V M = v£ ) +v M (y,t,t) (24) 

We switch on perturbations only in the z and x-directions (transverse directions): 

a>z(y, t, r), v z (y, t, r), a x (y, t, r), v x (y, t, r) 



5 



Since the Chern-Simons contribution to the current depends only on the intrinsic gauge 
fields on the cutoff surface, its flow is trivial. The non-trivial part of a possible flow is 
completely contained in the covariant currents 



j 



(l)x 



9 2 v 



9a 



-9* A 



(25) 
(26) 



r=A 



We define our transport coefficients as the response to the perturbations and in terms 
of the previously defined covariant currents as 



J [1)x = a C MEe{rtxyz)F$ v + a axial e{rtxyz)F^ A 



ir- 



J: 



(1)x = a axml e(rtxyz)Fff v + a 55 e(rtxyz)F^ A 



(27) 
(28) 



o 'axial defines the vector current generated by an external axial magnetic field. Observe 
that, in order not to have F}^y\ = ^ identically, one has to turn on the perturbations 
a x (y,t,r),v x (y,t,r). However, these do not play a role when studying the flow of the 
anomalous conductivities for they induce contributions that tend to zero in the low u>, 
low k limit, very much as occurs in 



The value of the background gauge fields is [21 



.4 



(o) 



oi — 



^r 2 H 



Vc 



(o) 



7 - 



fir 2 H 



(29) 
(30) 



The integration constants a and 7 can be fixed by e.g. demanding that the gauge fields 
vanish on the horizon. In any case the covariant currents do not depend on these integration 
constants. The consistent currents ( 1T9|) . ( 120]) do however depend on them through the 
Chern-Simons currents. For a discussion of this dependence see [20| [21] . 

The procedure consists of using the equations of motion to find the value of d r a, where 
a a generic conductivity defined at some hyper surface r = A. In fact, we only need the 
equations of motion projected onto x and the Bianchi identity associated with the indices 
(r, y, z) to obtain an expression of the derivative with respect to r of the different transport 
coefficients. From the simple form of our metric it can be seen that the vector normal to a 
hyper surface of x=constant reads = ^g xx (0, 0, 1, 0, 0). Hence, for the vector gauge field 
we have 



M 



9 2 v 



jtj jpNM ,0 MNPQRfpA pV \ 



(31) 
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Taking advantage of the relation VnF nm 
the currents (I25D and (1261). we arrive at 



7=^iv (y/~9^ NM ) an d the definition of 



(32) 



where we have neglected Fy^ tx ; Fy^ yx for these modes lead to vanishing contributions in 
the low momentum and low frequency limit, as mentioned before. Besides, we have carried 
out the contraction e xNPQR F$ P F% R = -4e rtxyz (F t ^F^ z + F^F? + F£F% + (Ao V)). The 
Bianchi identity to first order associated with indices (r, y, z) reads 



drF {vA}d) + dyF {v,A}(D + dzF {vAm = o 



Assuming d z Fi 



J 9rr9yyd z J y ^ = we obtain 



-g 



9{v,A} a v J {v,A} 



Now, making use of these ingredients, the computation of d r a is immediate: 



(33) 



(34) 



d r (JcME = hm 

ui,k— >0 



JX{1) 



e[rtxyz)Fy^ 



[e(rtxy z) F^ 1 ^ 2 



d r FYS 1] 



yz 



yz v c V ^y^)^- yz 
Plugging f )32|) and f )34|) into fl35|) we find, in momentum space 



a M =0 



(35) 



d r (J C ME = hm 



1 o ci-A(O) i ' ; 9zz9rr 2 



<? " e(rtxyz)Fy Z ^ 



Taking the limit cu, A; — >• and substituting F t 



V(0) 



— d r A Q 



(0) 



(36) 

2^T a -, we get the 



following flow equation for the chiral magnetic conductivity 



d r (TCME = -24k 



(37) 



whose solution is 



2 

&CME = C+ llK^f . (38) 

C is an integration constant that we must fix. In order to do that, we impose in- 
falling boundary conditions for the perturbations (or, equivalently, regularity at the horizon 
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r = Tn [EE])- This in turn implies that the fields must depend only on the combination 
dv = dt + ^J^dr [13]. Therefore, in the A r = gauge, we have 

d r A x = x [^d t A x atr = r H (39) 
V 9u 

This condition forces directly J^ x {r = th) to be 

jW^r = r H ) ~ F Wxr (r = r H ) ~ E i (40) 

Imposing these infalling boundary conditions results therefore in a vanishing chiral mag- 
netic conductivity at the horizon for the covariant current Thus the integration constant 
C can be fixed simply by the condition 

a C M E (r = r H ) = Q^C= -12k/i 5 = (41) 



(Recall that n 



24vr 2 / ' 



The transport coefficients ( 12"T|) and ( I2"gj) can be calculated in an analogous way. 
For the axial current the projected equation is 



which implies 



V7 rpNM i MNPQRt rpA t?A , j?V t?V \ 

— V N t A + — £ {£ NP t QR + * NP^QR) 

9a a 



(42) 



dr£ )x = y/^gUnT-v (^F(r + ^ 1)V ) (43) 
The values of the conductivities at r = A then read 



acME(A) = ^ ( 1 - £ ) (44) 



2tt 2 V A 2 

<W(A) = ^ (l " J ) (45) 



.2 



■""Note that the consistent currents might have non- vanishing chiral magnetic conductivity on the horizon 
due to the Chern-Simons contribution and depending on the value of the integration constant a. 
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As expected, the result in the limit r — > oo is precisely the one obtained in [21] using 
AdS/CFT techniques. 

In view of the topological nature and the non-renormalization theorem for the chiral 
magnetic conductivity it is at first sight somewhat surprising to find a non-trivial flow. This 
flow becomes however natural if we define the chemical potential in its elementary way as 
the energy needed to introduce one unit of charge into the ensemble. In the holographic 
dual this corresponds to bring a unit of charge from the boundary, now situated at r = A 
behind the horizon. The energy difference between a unit of charge at the boundary and 
a unit of charge at the horizon is just given by A (A) — Aq(th) = MA). This defines an 
effective chemical potential in the theory equipped with the cutoff A. In fact the definition of 
such an effective chemical potential is natural even in field theory. If we have a momentum 
cutoff of order A we can localize a unit charge only inside a volume within a radius of order 
1/A. Thermalizing this unit of charge means spreading it out over the entire ensemble. The 
difference in energy between the two configurations, the unit of charge localized within 1/A 
and spread out over the ensemble again is the effective chemical potential. 

All the anomalous conductivities can therefore be expressed in the form 

*A)-^A). (47) 

The are linear in the chemical potential and the numerical coefficient is independent of the 
cutoff. In this sense they obey the expected non-renormalization theorem. 

3.2 Inclusion of metric perturbations 

In this section we compute the flow equations for the Green's functions associated with 
generic response. The method can be described as follows: we need to consider two equations. 
One is the constitutive equation 



N 

\Oj) (48) 

i 

and the other one is the covariant holographic definition of the one-point functions, eval- 
uated on some perturbed state. Generically, these would be a functional of the perturbations 
and its derivatives (the dot means d/dr. We will be using both notations indistinctly). 

N 

w ■ n 'i°^ ( 49 ) 

i 

Taking the r-derivative in both equations, we can force them to be equal. Observe that, 
from ()49p . we expect terms containing H*^. After using the equations of motion, we will 
be left with some expression involving only (j) and <f>. Then, by equating (1481) and (1491) . it 
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is possible to find a formula for (fij = £\ K]4>i so that eventually we are able to write the 
r-derivative of as an expansion in the perturbations only. 

On the other hand, differentiating (EES]) and using again <pj = ^\ K^fa, we are lead to an 

expression in terms of Gj, Gj and 0j. 

Imposing that the r-derivative of ( 1491) and that of ( 148]) are identical, we finally arrive at 

N 

= ^]^ (50) 

i 

where A 1 - is a functional of G* and their first derivatives. Assuming now that the different 
perturbations are independent from each other, we get N independent equations 



A) = (51) 

which are nothing but differential equations for G*. Remarkable enough, the flow equa- 
tions for the retarded correlators are of first order in r-derivatives. 

3.2.1 Application to the anomalous conductivities 

In what follows we will derive the flow equations in the presence of a pure gauge Chern- 
Simons term (no gravitational anomaly) by using the procedure detailed in the previous 
section. The model reads 



S = S EH + S gh + Y^c[ <a ^ (~\f mn F mn + ^ NPQR A M F NP F Q ^\ (52) 

where Seh denotes the Einstein-Hilbert action with negative cosmological constant and 
Sqh is the Gibbons- Hawking term on the boundary r = A. The Chern-Simons coupling is 
here related to the anomaly for a single chiral fermion by k = —G/(2tt). 

Since we need now the precise equations of motion for the metric fluctuations we will 
specialize the analysis to a Reissner-Nordstrom AdS Black Brane 



ds 2 = T — (-f(r)dt 2 + dx 2 ) + ^rr^dr 2 (53) 
L z r j[t) 

A<® = <f>{r)dt = -^f (54) 

The integration constant in the gauge field is set that it vanishes for r — > oo. The horizon 
of the black hole is located at r = rfj and the blackening factor is f(r) = 1 — + ^ . The 

parameters M, Q are related to the chemical potential at infinity /i and r# by M = j^ + ^ir, 

2 ^ 

Q = ^S-. Finally, the Hawking temperature is given by 
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T 



4vrL 2 ' 



27rr 5 H 



(55) 



In what follows, we consider perturbations of momentum k in the y-direction at zero 
frequency. It is only necessary to turn on the shear sector, that is, the perturbations are 

x, z 0. It is more convenient to work with the coordinate 



written as a a , /i", where a 

r 2 

u = instead of r. 



The equations of motion for the perturbations derived from 
O(k), read 



when uj = and to 



where R = 



fW f(u) 

f^\ u )-!^M-3auB'(u) 
u 



ike a/ 3K 



Bp{u) 
/(«) ' 



(56) 
(57) 



The operators that we will be working with have the following form when evaluated in a 
perturbed state (for further details see [9]) 



5J a 
5t' 



H 



8itGL 3 
r 4 H f(u 



K - -K 



u 



8irGL 5 u 

where the prime stands for d/du. Differentiating ({55]) and fl59|) we are left with 



(58) 
(59) 



(6 J 
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8nGL 3 



(aMf(u) + a> a (u)f(u)-vtin 



8ttGL 5 



a 



m 
/(«) 



6 3f'(u) 



u- 



uf(u) 



(60) 
(61) 



In order to handle the terms, we evaluate the above expressions on-shell, yielding 



(SJ 



11 



8nGL 3 



8nGL 5 u 



m 
/(«) 



6 3/'(u) 



u- 



uf(u) 



Kike a pap 



3au , 
a Q 



(62) 
(63) 



1 At zero frequency the fields hy decouple from the system and thus will not be considered (see [8]). 
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Now, observe that, since h' t a = ^fj^St? + \h% and a' a = (^^-8J a + /ihfj jfa, §5§ 



turns into 



8nGL 5 u 



au 



/(«) 



u 1 



6t? + ^ 6 J a 



(64) 



Plugging the constitutive relations (e xz = 1) 



SJ. 



const 



G xx 5 a ^a p + G*V% + P xt ^h^ + P zt ^h% 
K const = G xx 5 a ?ap + G* 2 e a % + Pf + Pf e^h^ 

into (I64I) . the remaining equation for (££")' involves only <pi and G*. 



(65) 
(66) 



On the other hand, we can take the u-derivative of (!65l) - (l66l) explicitly and then make 
use of (j58|) - (l59p to end up having an equation in terms of 0j, G* and GJ. 

Finally, imposing (SJ" onst )' = (5J a )' and (St" const )' = (Stf)' and assuming that the 
perturbations 0, are independent from each other, we find 



G 



1 {{G xx f - (G xz f) - 8 f? 2 f^ (P xt G xx - P zt G xz ) = 



/xz , 167rGL 3 ^ TT ^ 8ttGL 5 u 

G 



f{u)r 



QXX QXZ 



H 



r 4 H f 2 (u) 



p ,xt + ( _^ + 



/(it)r 



2 



(p*'Gf + i> 'a 1 /') = 

8ttGL 3 
f(u)r 2 H 



~8nGL 3 



hzik 



8ttGL 5 



+P 



xt 



p,zt + G xz(_ 



/ 8vrGL 5 
8ttGL 3 



r A H f 2 (u 



pzt \ pzt 



U 



pxt _ 



f{u) + 3* 
u 



px 



xt 



/(«) 

8ttGL 5 u 







+P 



8ttGL 5 



/ 2 («) 



pxt pzt _|_ 



3 />) 8nGL 3 
P{u)t% e + u /(«) + /(«),* 



G a 







(67) 
(68) 



(69) 



(70) 
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riixx i 8irGL 3 lnxxnxx n xzrixz\ 8irGL 5 u xtnxx 8irGL 5 u ztnxz _ 

= _ Gr Y3 +/uG - (71) 

V M /(")/ 

r~ilxz I 8ttGL i (~ixx s~ixz , /~<xz /~ixx\ 8ttGL U / ryxtr^xz rjzt /~ixx\ 



/(u)rj, ; /(ajri < ' V r'„P(u) ' u f(u) 

r HJ \ u ) \ u J\ u ) J 8irGL b u \ u l ) 



J(u)r 2 H ^ ) f(u)r 2 H 6 / 2 Hr|, 



By directly studying the structure of the solutions to (j56"]) - (j57]) . it can be realized that 
qxx _ pxt _ qxx _ q f or w _ g an j ^ g rs ^. orc [ er j n fa Furthermore, all the anomalous 

correlators are of order k or higher. A more detailed study of (|67p - (l74p is left for section 

3.3 Flow of the transport coefficients as two point functions 

As suggested in Section 2, we could have determined the flow by simply considering the 
system to be restricted to live between the horizon and a cutoff surface placed at A. It is 
hence expected that the transport coefficients at the boundary can be computed by finding 
the corresponding 2-point functions. The boundary value of the perturbations, whose bulk- 
to-boundary propagator is normalized at the cutoff, work as the sources for the different 
operators of the dual theory. 

Henceforth, the perturbations will be rearranged in a vector It is more conve- 

nient to use the Fourier transformed quantity 



*(«,*")= / 7^uHHe- M * S (75) 



d d k w 
(2tt) £ 

The explicit expression for is 

*J(u) = (B x (u),h x t (u),B z (u),h z t (u)) (76) 
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being B a = a a /fi. To proceed, one can follow [22] and assume the general form of a 
boundary action 



5S^ 



r=A 



d d k 

(27r)° 



(77) 



In order to get the solution of the system (I56 i) -( l57l) to first order in momentum we expand 
the fields in the (dimensionless) quantity p = Hence 



u 



h^ a (u)+ P h: 

= B£\u)+pB i 



,u 



(l).a. 



Mi 



U 



U 



(78) 
(79) 



The system can be solved perturbatively . To calculate the retarded correlators at r = A 
(or, equivalently, at u = u c = r 2 H /K 2 ) we only need to solve the equations for the per- 
turbations with infalling boundary conditions, on the one hand, and boundary conditions 
$^(« c ) = (ftl on the other [9]. This procedure should give us the desired Green's functions, 
after taking the variation of (ITT)) with respect to the fields at u = u c (which act as sources for 
their corresponding operators). Recall that, as explained in Section 2, the bulk-to-boundary 
propagator must be normalized at r = A, that is, if we have 



$i(«) = Fj(fc, M )0 



then Fj(k,u c ) = 1. Notice that the relation between the boundary value at u = u c and 



Fj(k,u c )(j)'l , so that the solution is preserved by these 



that at u = is simply <$ k ^ 
manipulations, as pointed out by [13] and [16]. The retarded two-point functions, from which 
we are able to read directly the transport coefficients, then have the form 



Gjj(k, u c ) = -2 lim (a im {Ff(k, u))' + B 
Where the Au and Bu matrices are [9] 



A 



r 



1 



— ^——Diaq I — 3af (u), — , — 3af (u) 



(82) 



B 



AdS+d 





"0 


3a 








4 

r H 





3 

V? 



















3a 













3 



(83) 
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Using again the the effective chemical potential 



MA) = li ( 1 - -§ ) , (84) 



the result for the anomalous correlators is 



/xt*at*\ W Kk fi t h\ -ikfi(A) 

{U 5J)= 2^G\ 1 -^) = (85) 

(5J*5tl) = {<H*6J*) = (l - ±y = ^! (86) 

^^-Mr^j - 12^ (87) 

Since limA->oo A*(A) = /i, these correlators coincide essentially with the ones derived in [9ffi 



3.3.1 Compatibility with the flow equations 

The system of first order differential equations fl67|) - f!74|) must be compatible with the result 
fl85|) - flH7|) encountered in the previous section. In order to check that it is so, the dissipative 
correlators play an important role. In the case tu = and to O(k), they read3 



£<XX pXt QXX Q 

SttGLV [) \f(u) u) m 

This solution implies that G xx and P xt = G xx are of order u or higher, whereas P xt 
contains a part which is of order O(k ,u°) (contact term). The remaining system, after 
substituting ( )88|) . (189|) and assuming that all the anomalous correlators are at least of 0(k), 
turns out to be (up to order k) 



3 The minus sign found in with respect to the result of [3] is due to the fact that in this reference the 
correlator that is studied is (5J a St t b ) 7 that differs from (<5J a (5^) by a factor of (b represents a spatial index) 
g tt g bb = -f(u) -> -1 at infinity. 

4 The limit P^iu = 0) is not well defined because we have not included the corresponding counterterms 
in ([55)1 . (|59|) . The reason is that they do not affect the anomalous correlators. 
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G' xx = (90) 

G "° = -fSv mk (91 > 

P' xt = (92) 

P' zt - fiG xz = (93) 

G' xx = (94) 

G ,xz = fiG xz (95) 

Pf = -Pf (-- f -M)- jr^-fM - (96) 

\u f(u) J 8tvGL 5 u \ u 2 J 

P' zt - nGf = [iP zt (97) 

Equation ( )96|) is in agreement with ( 1891) . In the end, the 2-point functions associated with 
dissipative transport coefficients decouple completely. Regarding the anomalous correlators, 
the above system of equations can be integrated easily, leading to 



GXZ = ^GT^ k{l - u ^ (98) 

2 

P Zt = GT = -Vj^JJ^k (1 - u c ) 2 (99) 

pf = ^^^(i-- c ) 3 (ioo) 



which is the same as (185l) - (187j) . The role played by the Chern-Simons term in (1521) is 
crucial to ensure that G xz presents a flow, for in its absence all the anomalous 2-point 
functions identically vanish. 



4 Gravitational Anomaly 

The study of the effect of the Gravitational Anomaly on the definition of the holographic 
operators is a non-trivial task, for the term A A R A R has not a well defined Dirichlet 
problem. This makes, strictly speaking, not possible to define generic operators. In [8], the 
problem was circumvented by arguing that any possible contribution vanishes asymptotically. 
However, now we are interested in the value of the transport coefficients at finite cutoff A, 
and therefore it is necessary to face this issue. 

4.1 The Model 

The four dimensional axial gravitational anomaly is induced holographically by a Chern- 
Simons term of the form [8] 
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Sacs = ^ J d 5 x^e MNP ® R A M R A BNP R B AQR (101) 

This action contributes to the boundary axial current as expected for a mixed anomaly. 
The complete action reads 



S = / d b x^g 

16nG 1 



R + 2K c - l -F MN F MN 



Sacs + &aem + Sg + Scsk (102) 



Where 



Saem = [ d 5 x^e MNP ® R A M F NP F QR , (103) 

4o7rG J 



s ° = ~l^^- hK < < 104 > 



t) 



Scsk = ~7Tn f d A x^hn M e MNPQR A N K PL D Q K L R . (105) 

ZnLr JdM 

Adding Scsk ensures that the anomalous Ward identity for gauge transformations de- 
pends only on the intrinsic curvature tensor on the boundary at r = A [8] . 



Indeed, the covariant current turns out to be 



16nGJ A = n B [F AB - 8e BACDE XK CF D D K B ] r=A (106) 
with a purely four dimensional divergence that on shell evaluates to 



D.J" = l —e° 



(107) 



167TG 1-3 ° P 9T (4)6op-"-(4)agr 

where e opqr = e nopgr is the four dimensional epsilon tensor. 
The bulk equations of motion are 

1 1 

Gain — ^ c Qmn = -FmlFn L — -F 2 gMN + ^^lpqr{m^b (F pl R b n ^ qr ) , (108) 
V N F NM = -e MNp Q R (KF NP F QR + \R A BNP R B AQR ) , (109) 

4.2 Contribution of the Gravitational Anomaly 

If we vary iSUcs> we are left with a term which spoils the variational problem 
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If we looked for a suitable counterterm to render the Dirichlet problem well-posed, we 
would end up finding Scsk- Indeed, this boundary contribution was firstly conceived as 
an analogue to the Gibbons- Hawking- York term. However, after varying S^cs + Scsk one 
realizes that the result 



is still problematic. Even worse, f illip can not be canceled easily, for, for instance, the 
ansatz 



is automatically zero. Thus in principle, there is not a straightforward way of having a 
well defined variational problem for this system. 

On the other hand, as aforementioned, we need Scsk, to have a four dimensional anoma- 
lous Ward identity at the boundary, so we will keep it. A hypothetical generic counterterm 
(if it exists) capable of solving all the problems, would probably ruin fllQ7p and therefore, by 
physical means, should not be considered. 

Even though the variational problem is not well-posed, we will still be able to derive the 
equations of motion by means of the analogue of the Euler-Lagrange equations for higher- 
derivative theories. The difficulty therefore reduces to the question How to treat 
holographically? Note that in [22J it was implicitly assumed that the Dirichlet problem is 
correctly defined, so we should go a little bit further in this case. 

Specializing for the shear sector, which is the one that interests us, and at second order 
in perturbations, fillip reads 



Other possible terms would vanish in the background fl53l . The strategy would be the 
following: Since f)113p does not affect two point functions involving only energy-momentum 
tensors or only currents, we know how to compute the correlators (T X T^) and (J X J Z ). (11131) 
only plays a role when calculating (TfJ z ), (J x Tf), and hence those are the ones for which 
the discussion of [22] does not apply. 

Following the method detailed in Section 3, it turns out that, taking only into account the 
gravitational anomaly 




(111) 




(112) 




(113) 
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(r,r) = o (ii4) 

(If 2?) = " i2 Mc)r2 (115) 

(note that we have directly substituted the value of A for a single left-handed fermion 
\/G = — 4g-)- The above results point again towards an effective — u c ). Therefore, by 
physical grounds, we expect the appearance of an effective temperature also. Note that the 
flows of the effective quantities must be consistent in the sense that they must be the same, 
no matter what correlator we are focusing on. Equations flll4p - flll5l) hint at the existence 
of an effective temperature for the system; this temperature does not flow with the cutoff 
scale, being always identical to the Hawking temperature. This conclusion is in agreement 
with the asymptotic values of [SJ. 

So we resolve that (11 111) must be treated in such a way that (TfJ z ), (J x T t z ), at finite 
cutoff, are consistent with a non-flowing temperature. 

It turns out that the method to achieve it is precisely the one that one would anticipate by 
general considerations: Taking advantage of the fact that the equations of motion for the 
shear sector 



h a> (u) 

hf(u) - -L±-L - 3auB'(u) + i\ke a g 
u 



(24au 3 - 6(1 - f(u))) 



B p {u) 



u 

+ (9au 3 -6(1- f(u)))B' p {u) + 2u(uhf (u))' , (116) 

fw /(«) 

3 - (2 , , _ „ 3 \ or _B p (u) 



+ik ^ \ " :) + 3 " ) h " + R iU) ' (117) 

happen to be of second order in derivatives (where A = ^—), we can solve completely 
the evolution as we did in Section 3.3 (imposing in-falling B.C. at the Horizon and Dirich- 
let B.C. at the boundary). Once the solutions are known (see the appendix), f illip will in 
general give a well determined surface contribution (when evaluated on-shell) that must be 
taken into account to calculate (TfJ z ), (J X T^). The result so obtained presents no flow in 
the temperature part. 

To be more concise, the boundary term fillip to be considered has the following form 



ik\r 2 H e afi 



<) 



2ttGL 4 

whose contribution, up to first order in k, is summarized 



uf(u)ap(k)ti t a (-k) (118) 
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ik\r 



2ttGL 4 



f'\u) (o) 
U— — rO 

d /(«, 



(*)^ 0) (-fc) 



(119) 



(Notice the factor ~ jk^ introduced to normalize the perturbation (see the appendix)). 
So the effect of (11 111) on the Green's functions can be reformulated as a modification, pre- 
scribed by (11191) . of the Bjj matrix. 



Even though (11191) only affects the correlator ^T^J^), S^cs + Scsk induces automatically 
a non-vanishing value for the components An 



A^ 2 of the matrix A. These contributions, 



which are perfectly treatable within the framework of [22], give rise to a correction of ^ J a Tf 

which is precisely of the same form of the one implemented by (I119p . As will be mentioned 
below, this turns out to be sufficient for the consistency condition (I127p to hold. 



The final form of the matrices Au and Bjj after implementing the shift driven by the 
Gravitational Anomaly is given by 



A 



ii 



16ttGL 5 



-3af(u) 






l 

u 



AiXkL 
' H 



uf(u) 



iSXkL^i 



U 






-3af(u) 




4iXkL 



uf{u) 



i&XkLu 
5 11 

r H 



1 
u 



(120) 



B 



AdS+d 



11 



IQitGL 5 



-4iXkL 







9au 3 -6(l-f(u)) 



3a 

3 






AiXkL 





9au s ~6(l~f(u)) 









3a 



'1211 



D, 



dCS 



H 



16ttGL 5 









4iXkL „ J'(u) 2 

n 11 










4i\kL u f'{u) 2 



' II 



/(«c) 



(122) 



The resulting anomalous 2-point functions are 
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ikK/j,(l — u c ) u(A) 
ik*. (1 -u c )V (-2 + a) 2 r| <V(1 - u c ) 2 T 2 



= ""»^_ = -ik^- (123) 



^Z%^ - * + 24 J < 124 > 

/7 ,«v (l-^c) 2 /i 2 zfcA(-2 + a) 2 rl ,/ ^ 2 (l - n c ) 2 T 2 \ 

^£-^)V v ikX (-2 + a) 2 r« _ //i 3 (l - Mc ) 3 , /i(A)T 2 

\ J t J t / — r [l- U c ) fi 777IZ ~~ ~ lK 



6Gtt v c/f ^ GL% V 12 vr 2 12 

(126) 

Observe that it is straightforward to verify that equations fll23l) - fll26p are compatible 
with the asymptotic value computed in j8]. Notice also that the temperature part re- 
mains constant as we move the boundary. The flow of the different correlators is con- 
sistent with respect to each other and the hypothesis of an effective chemical potential 

/i(A) = )jl (l — = /i(l — u c ) is reinforced by the results extracted from the terms propor- 
tional to A. 



5 Discussion and Conclusion 

We have studied the holographic cutoff flow of the anomalous transport coefficients. This 
has been done by defining a bottom up model that implements both the axial and the mixed 
axial-gravitational anomalies. The flow has been studied by analyzing the dependence of the 
anomalous Green's functions on the radial position, A, of the boundary. We have presented 
several prescriptions to compute such flow and finally obtained it by adapting the method 
implemented in j8], [9] for the case A — > oo. 

It is a remarkable fact that the chiral magnetic conductivity suffers from a flow even in the 
non-backreacted case. In fact, this could have been anticipated by noticing that regularity at 
the horizon imposes that in the deep IR the constitutive relations are only compatible with 
an electric conductivity ([13]), so that if a system exhibits a chiral magnetic conductivity in 
the UV it must be due to a non-trivial flow. 

When considering the gravitational anomaly, a Dirichlet boundary condition is not enough 
anymore to define the variational problem properly. A generic definition of suitable opera- 
tors, if any, therefore requires further discussion in this case. In this paper we have simply 
focused on computing 2-point functions, without discussing general definitions of the corre- 
sponding operators. The term which spoils the variational principle has been dealt with by 
considering its effect on the on-shell action. This procedure, which can be seen to be the most 
natural one by using physical arguments, yields 2-point functions that are consistent and 
whose flows do not get in contradiction with the result found in the absence of gravitational 
anomaly. Moreover, in the spirit of [22], that the matrix of correlators Qu obeys 
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(127) 



represents a non-trivial consistency check. 

The result f ll23l) -( |T26l) shows that the temperature remains constant (Hawking temper- 
ature) whereas the chemical potential presents a flow that is easily interpretable in terms of 
the energy necessary to bring a unit of charge from the horizon to the boundary. Observe, 
however, that all the correlators are written for a metric with g tt ~ — r 2 /(r), and hence there 
is an implicit redshift factor between observers living in one hypersurface placed at r = A 
and another one at r = A'. 

From the point of view of the boundary theory, these outcomes indicate that the pure gauge 
Chern-Simons term does not affect the boundary operators but influences the anomalous 
correlators through the flow equations, forcing them to have a non-vanishing value at the 
boundary, whereas the gravitational-gauge Chern-Simons term happens not to have any im- 
pact by means of the evolution equations, but to induce new covariant contributions, that 
are first order in k, to the operators, so that the constant T 2 part is present at any value of 
the r-coordinate. 



Appendix: 

Solutions at zero frequency and normalized at finite cut- 
off u c 

Case A = 



B a (u) = B a + H a (u-u c 



2(1 + 4a) 2 (-l + u c {-\ + au c )) 
x((l + 4a) (u - u c )(H> 3 + H p u c + a{?>B p {2 + u c ) + H P (A - u c (2 + 3u c )))) 



x 



-2y/l + 4a(-2 + a(-2 + 3u))(B fi - H p u e )(-1 



u c (— 1 + au c ))(ArcTanh 
+ArcTanh 



-1 + 


2au 




\- 4a _ 


"1 - 


2au c 


Ivi 


+ 4a. 



)) 

(128) 
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XI 



2(-l - 4a) 3 / 2 (-l + 0(-l + «c(-l + a« c )) 2 
-2(-l - 4a) 3/2 iJ Q (-l + u(-l + au))(-l + u c (-l + au c )) + 

+khe aP {-i (y/-l - 4a - iVl + 4a) ^(1 + it)(l + w c ) + 



+ 



+a 2 3B^ ^2iV-l - 4au 2 + zV-1 - 4auu 2 + Vl + 4au 2 (2 + u c ) 
2 if' 3 ^2zV-l -4a(2 - 3u c )u 2 + zV-1 - 4au(4 - 3u c )u 2 + Vl + 4au 2 (4 - u c (2 + 3u c )) 

-3iaB l3 (2V-l - 4a - 2iy/\ + 4a + 2V-1 - 4aw c - iy/l + 4au c )- 
-ZiauB 13 W—i - 4a - 2iy/l + 4a + - 4att c - iVl + 4aw c j + 
+aH ,3 (-AW-l - 4a - 4Vl + 4a + Vl + 4au 2 (l + « c ))+ 



2zV-l - 4a + 2Vl + 4a + 7iV-l - 4aw c + 3 Vl + 4aw c J + 

+uaH ,3 {-4iV-l - 4a - 4Vl + 4a) + 
+uaH l3 u c - 4a + 2Vl + 4a + 4zV-l - 4au c + 3Vl + 4au c )) + 

+6mfc(— 1 + it(— 1 + au))Re a p(B^ — H^u c )(—1 + u c (—l + aw c ))ArcTan 



2au 



+Qak(— 1 + tt( — 1 + au))K€ a p{BP — H^u c )(—1 + w c (— 1 + aw c ))ArcTanh 



V-l -4a 
1 — 2aw c 

Vl + 4a 
(129) 



+ 
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Case k = 



B a (u) = B a + H a (u - u c )+ 



-2 + a)ku c e a p(- 



2i(-2 + a(-2 + 3w))ArcTanh 



\—2au 



X 



a 



6(2 -a)a 3V ' ' -'""""^v (1 + 4a) 3 / 2 

x + a(3(l + a(7 + 2a(7 + a)))B^ + 4(8 + a(2 + a) (9 + 2a))H p - 

-3(1 + a(7 + 2a(7 + a)))H p u c )) + 

-2aVl + 4a(u - n c )(6a(i? /3 - tf^-8 + u)) + 6a(i^ - ^(-8 + u))u c - 8aH^u 2 c - 



-3a 2 (3i? /3 (-4 + «)(1 + w c ) + (u(10 + 3u) + w(10 + 3u)u c - 2(-8 + u)u 2 - 4(7 + 6u c ))) + 

+8# /J (l + u c )+ 

+a 4 (^(12 + u c (18 + (-59 + 12u(2 + 3u))u c )) + 3^(2 + w c (5 + 12(— 1 + u - u c )u c )))+ 

+a 3 9B> 5 (4 + u(-4 + (-4 + w c )w c ) - w c (-5 + w c + w 2 )) + 
+a 3 F /3 (29 + (23 - 72u c )u c + 9u 2 (-4 + (-4 + u c )u c ) + 6u(-4 + w c (-4 + 5u c )))) + 
+ (-2 + a(-2 + 3u))(-l + u c (-l + au c )){2(AH p + a(3(l + a(7 + 2a(7 + a)))B p + 



1 — 2aw c 



+ 



+4(8 + a(2 + a) (9 + 2a)) H p - 3(1 + a(7 + 2a(7 + a)))#V))ArcTanh 

. VI + 4a 

+ (1 + a)(l + 4a) 3/2 (-4ff /3 + a(-3i? /3 - AH P +3H^u c ))x 

x(Log[-l + u(-l + au)} - Log[-l + u c (-l + au c )])))) 

(130) 



where ° _ ( _ 1 _ 4a)3 /2 1 (1+Uc _ au 2) 
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H?{u) = H< 



k (1 



U 



aw 2 ) u c x 



(-1 + u)(-l + u(-l + au)) (l-tt)e 

(-l + u c )(-l + u c (-l + au c )) + 2(-l -4a) 3 / 2 a 2 

1 

x (— - — — — — — — 2aVl + 4a(u — u c ) x 

v (-l + w(-l + au))(-l + « c (-l + au c )) 

x (4H P (1 + u)(l + m c ) + aB p {3 + 5u + 5(1 + u)u c ) + 

+aH 13 (18 + w(25 + u) + 22w c + w(25 + u)u c - 4(1 + u)u 2 c ) - 

-3a 2 i^ 3 (-5 + 2w 2 (l + u c ) + u c (-7 + 2u c ) + u{-7 + 2(-3 + OO) + 



ifV (18 + 24u c - 22u 2 - 6u 3 (l 



u(39 + 5(7 — 5u c )u c 



u 2 (1 



u 



;)) 



+3a 3 I? /3 (1 - 8u 2 - uu c (3 + 8u c ) + 2u 2 (-4 + (-4 + u c )u c )) + 
+i^a 3 (4 - 4u c (2 + 5u c ) + u 2 (-20 + (-20 + u c )u c ) + 6u 3 (-4 + (-4 + u c )u c ))- 
-H^a 3 u(5 + 7u c (l + 5u c )) + #V (4m + Au c + 6ww c + (-2 + u(-5 + 4w(5 + 6w)))w 2 ) + 

+£V(2u c + u(2 + m c (5 + 24ww c )))) + 
+2(4^ + a(3(l + a(7 + 2a(7 + a)))B p + 4(8 + a(2 + a) (9 + 2a))H p - 



-3(1 + a(7 + 2a(7 + a)))H p u c )) I ArcTanh 



2au 



ArcTanh 



2aw r 



VI + 4a J La/1 + 4a. 

+ (1 + a)(l + 4a) 3/2 (-4# /3 + a(-3fi /3 + i^(-4 + 3u c ))) x 

x (Log[-l + u(-l + au)} - Log[-l + u c (-l + au c )])) 

(131) 
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